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General 


The  orientation  distribuiton  function  F  of  nonspherical  aerosol 
particles  is  very  decisive  for  many  cloud  properties  such  as  light  scatte¬ 
ring  and  radiative  transfer  (1,2),  diffusional  transport  mechanics  (3-5) 
and  average  rate  of  sedimentation  (6). 

It  is  also  very  important  in  rheology  for  the  constitutive  properties 
of  non-Newtonian  fluids  (7-9). 

In  general,  the  orientation  function  is  affected  by  two  opposing  physical 
factors:  The  randomizing  action  of  the  rotational  Brownian  motion  (or  micro- 
turbulent  eddies)  and  the  orienting  influence  of  flow  gradient. 

Up  to  now,  the  studied  situations  were  those  of  the  (asymptotic)  weak  and 
strong  field  in  which  the  rotational  Peclet  number  <*  which  indicates  the 
rati  o  between  the  above  factors ,  vi  z .  =  °/  x>  ^  was 

1  or  ot  i  .  Likewise, 

the  investigations  were  rather  concerned  with  the  simple  shear  flows. 

Aim  of  Study 

Since  many  cases  of  significance  are  characterized  by  values  of  cx 
of  the  order  of  unity,  as  occurs  in  the  free  atmosphere  and  with  typically- 
sized  particles,  it  became  of  interest  to  study  these  cases.  In  addition, 
the  general  laminar  field  where  the  nine  components  of  the  gradient  tensor 
had  to  be  treated  once  for  its  own  significance  and  second  time  for  employing 
it  in  models  for  the  real,  turbulent  atmosphere. 

Basic  equations  and  assumptions 

The  basic  equation  in  the>  reported  study  was  the  (source-free)  Fokker-Planck 
equation  of  conservation  in  angle-space 


c)F/bt  +  ?•(  Fu  -  [1] 


For  nomenclature  see  Appendix 


-  2  - 


We  have  considered  a  system  of  small  spheroidal  particles  immersed  in  a 
general  laminar  flow  which  is  given  in  their  vicinity  by  the  linear  relationship 


U  (*,t  )  -  W  (t)  •  x  .  [2] 


The  small  particles  simulate  (asymptotically)  straight  long  fibers  or  flat 
platelets  which  are  used  in  many  applications. 

The  rotational  velocity  of  the  particles  cj  of  Eq.[l]  was  taken  to  be 
given  by  the  Jeffery’s  equation  (10) 


Mi  =  iHvVj +  s'‘i  O]  im 

where  l  }  ^  ,  k  ->  1,2,3  are  cyclic  permutation  indices,  a,  ,  «■,  ak 
are  the  semi- axes  of  the  el  1  ipsoids ,  "V  '  ~  W'-  W.'  *£'.=  W*  ■+■  w'  . 

M  y  y  y  * 

=  3  u‘  /  dx'j  and  u’  x  are  respectively  the 

(fluid)  velocities  and  location  vector  in  the  body-locked  coordinate  system  x' 
(Fig.l ) 


1 
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Here,  using  the  usual  kinematic  relationships  between  0  ,  ^  =  q  ) 
and  co  '  co'  to'  ,  and  the  similarity  transformation  between  the  *’ 
and  x  systems,  we  finally  obtained  that 


*  =  V/2  +  ctg  (  £i{  cf  ~£Z[  S'?) 


-  ^  (  ^32  c;  f  -  Q.(s2'i)) 


and 


[4] 


6  ^ii3i  c  6  &  mi  s  o  +  -4-  s  4e  ( (3. 


+  +S„s4*) 


in  which  Eib  J.  (-v-  <  _ 

V  VlK  +  A  )  >  (?  _  c  ^  +  ^ 


kf  i  * 


(A  ,  9  ( Wi2  -  W,3  j  >  n  -  J_  (  u/ 

Q=  2<w  „  +w„-2Wll)  VIK  =  Wi.-W 

J  ki  , 


Wck“  W*.  ,  Wt 


h  ^U,  /  )Xk  >  -  COS  -m  -f  }  c^f  =•  c.os'^v 


•  n  -vn  -f  ,  -f  = 


'o'*™*  etc.  ,  X  -  (  R2-  0/(Ra-M)  ,  R  =  Q,/a; 


Now,  after  non-dimensionalization  by  designating  $ik  =  W;^/  W0 
and  I  -  t  W0  (  w0  -  a  typical  gradient  component),  we  could  convert 
the  Fokker-Planck  equation  to  the  form* 


; - - - 

The  sitfr,  is  dropped  out  for  simplify  sake. 
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ar/tt  -  ,,"<1  +  ~a*  <F^  =  a  F/“ 


(5) 


where  D  is  the  mid-diameter  rotational  diffusion  coefficient  of  the  axially  symme¬ 
tric  particles.  The  normalization  condition  was  taken,  as  ususal ,  to  be 


2ir  it 


C6] 


O  O 


Solution  of  the  Fokker-Planck  Equation 

According  to  the  theory  of  Fourier  series,  we  wrote  down  the  solution  to 
Eq.  [5]  as  a  series  of  spherical  harmonics 


V 

/. . 

in=o  y*-o 


l  o-  f *  j r*o+)c  * 


w, 

f\  (t) 

irt 


Thus,  inserting  [7]  into  [5]  and  taking  account  of  the  symmetry  proper¬ 
ties  of  the  physical  problem,  we  obtained  for  the  A-coe^ficients  the  set  of 
the  simultaneous  equations 


d  _ 


-n+a. 


ol  + 


«n 


L_  >_o  -  ) 

j=n-i  f=»--2  (-1 


x  [iy*  fT  +  c-  ok  c,r*  a  f  i 

(,r  f>r  '-k^J 


[3] 


prr) 

^  > 

and  the  size  parameter  x 


prri 
) 

defined  above. 


in  which  .  C  (  ,  j -*  are  functions  of  ,  *  >  j  a,jQl 


The  latter  were  numerically  solved  by  the  employment  of  a  fast  differential 
equations-solver  based  on  an  extrapolation  method  of  Burlish  and  Stoer.(ll) 

It  turned  out  that  our  general  solution  included  previous  work  as  special 
cases. 

*See  Proceedings  of  the  1983  CSL  Scientific  Conference,  Aberdeen  Proving  Ground , 

Md.  (to  be  published)  • 
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Numerical  Results 


Considering  realistic  situations  of  particles  with  a  characteristic 
size  of  r  io'f-  I  o’5  crn.  (whose  computed  D  is 

£>  ^  j  hc'1  )»  typical  flow  velocities  and  gradients  of  1-10 

m.sec"^  and  1-10  sec~\  respectively,  we  came  out  with  Peclet  number 
values  of  c*.  ^  1-10. 

As  indicative  examples,  we  bring  out  here  the  cases  of: 
i.  A  point  source  flow  whose  velocity  is  given  by 


a  -  ^0  r/  w  r3 


and  whose  gradient  is 


3  X,'  y  • 

r«r  * 


[9] 


noj 


where  x;  is  normalized  by  a  typical  length  r„  t  r  (x  ) 

is  the  radius  vector,  v/e  =  ’\»/Jiivr03  and  a  constant. 

ii.  A  round  laminar  jet  whose  axial  C  u. ,  ")  and  radial  (ur) 
velocity  components  are  (12) 


u  , 


*  l 


[ll] 


and 


-  iV'OjO) 


[12] 


and  the  gradient  components  are 


-  6  - 


Wh=  -  q  O-3iV<0,  w lr=r  -  fq  I  f  =  W„  |/)r 


W 


-  =  -f  *i  ( '  -  1  r  + 


)  -i 


Ur 


[13] 


2  l  v  *  •*  »6 

where  i  =  [f/X,  ,  {.($)-  (  '  ■+-  ^V^)-2,  y  is  determined  by 


the  jets'  momentum  J  ,  7 


(16/3)  TT  | 


'*  **  .  r 

distance  (again)  and  r  ,  x;  are  later  on  non-dimensional ized  by  r 


is  the  radial 


Indicative  results  for  these  two  significant  flows  are  presented  in  Figs.  2,3 
and  4,5  in  which  a  preferred  orientation  of  the  considered  particles  is  clearly 
seen  by  the  pronounced  maxima  of  F. 


8 


Fig.  2.  The  orientation  dsitribution  function  F  vs.  the  (Culer) 

angle  9  in  a  point  source  flow.  «  =  1  ,  r  -  0.02 

(platelets)  and  R  =  50  (straight  fibers). 

t  (  -  t  W, )  ^  1 }  x,=  A-  +  =o  (random  orientation)  ■ 

1  R  -  S  O'  f  =  O ,  z  4  --  n/M  ;  —  —  ,  R  =  o.Ox :  3  -A-0  ^  -4 


ttA  . 


Fig.  3.  A  perspective  view  of  particle  orientation  in  a  point  source 

field  (schematical ) . 


Fig.  4.  The  orientation  distribution  function  F  vs.  the(Euler)  angle  e 

flow  of  a  laminar,  axi -symmetric  jet.  I  ,  f?  =  o.oi,  0.2  ,5,  100 

/\> 

f  =  ( t  Wo)  z  ,  xi  ~  1  i  *z  =  *3  )  i  =  0.  s  ,  J  -  i  ,  f  —  oj  rr/i| 

'  R=  °°'  >  *  °>  2  -  0.0»jf  =  n/i,  .3  -  R  =  0.a  jts0. 

TT/A  •  S-R  =  S,  -f  -  rt/^  ^ 

7  -  K  •.  100.  f  0  j  8  -  *  ---  fOO,  -f  =.  TT/^  , 


-  8  - 


X  3 


xa 


X, 

Fig. 5.  A  perspective  view  of  particle  orientation  in  a  laminar  jet  (schematical ) 
Summary 

The  orientation  distribution  function  of  fibrous  and  platelet-like  aerosol 
particles  was  calculated  for  a  general  laminar  flow  by  solving  the  Fokker-Planck 
equation. 

It  has  been  found  that  this  function  shows  maxima  and  minima,  which  indi¬ 
cates  a  preferred  orientation.  The  cases  of  a  point  source  and  a  lani,.ar  jet  are 
brought  out  as  an  example. 
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Future  Plans 


Based  on  the  reported  study, the  next  stage  of  research  on  the  orientation  within 
a  turbulent  cloud  is  now  in  progress. 
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Appendix:  Nomenclature 


Qi  }  otj  f  ak  -  half  axes  of  an  ellipsoidal  particle 

,  f)  -  coefficients  in  expansion  as  spherical  harmonics 
-  coefficients  in  Eq.  [8] 

D,  D  -  rotational  diffusion  tensor  and  coefficient,  respectively 

defined  in  text 

F  -  orientation  distribuiton  function 

1  -  indices 

r  -radius  vector;  r0  -  normalizing  radius  vector 

R  =  a,  /  -u 

5 it,  ,  -  defined  in  text 

t  -  time 

u  -  flow  velocity 

v'lk  ,  -  defined  in  text 

wlk  >  W  ik-  -  component  of  the  gradient  tensor;  -normalizing  component 
x',  x  -  location  vectors. 


Greek  Letters 

-  rotational  Peclet  numbers 

^  -  parameter  related  to  jets  momentum 

-  Kroneckers 1  delta 

-  Eulers'  angles 


o>  -  rotational  velocity  of  the  particle 
v  -nabla  operator  in  angle-space 
A  -  Laplacian  operator  in  angle-space 


